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ABSTRACT 

Relativistic integral representation in terms of experimental neutron- 
proton scattering phase shifts alone is used to compute the charge form 
factor of the deuteron Gcd{Q^)- The results of numerical calculations of 
\Gcd{Q'^)\ are presented in the interval of the four-momentum transfers 
squared < < 35 fm^^. Zero and the prominent secondary maximum 
in \Gcd{Q'^)\ are the direct consequences of the change of sign in the 
experimental ^Si- phase shifts. Till the point ~ 20 fm^^ the total 
relativistic correction to |Gc£;(Q^)| is positive and reaches the maximal 
value of 25% at ~ Ufm-'^. 







Deuteron is the brightest example of intersection of nuclear and particle physics. During 
more then sixty years it serves as source of important information about the nuclear forces, 
mesonic and baryonic degrees of freedoms in nuclei, rclativistic effects and a possible role of 
quarks in nuclear structure. Therefore it is not surprising that currently the electromagnetic 
(EM) structure of the deuteron is a subject of intensive theoretical (the list of publication 
is immense) and experimental investigations. 

With new experimental data from Jefferson Lab on elastic electron-deuteron scattering 
expected in the near future §], at momentum transfers in the GeV-range, one needs 
to develop relativistic approaches to the (np)-bound state problem. Recent experimental 
results from MIT-Bates |^ provided the first experimental evidence for a zero in the deuteron 
charge form factor Gcd at about Q^= 20 fm~^ predicted in a number of theoretical models 
(or not predicted, as in some kinds of quark models). Here we report new results of numerical 
calculations of Gcd ■ These calculations are based on the approach to the relativistic impulse 
approximation, which was briefly discussed in ref. Q (see also the review and, especially, 
the references herein). The more detailed formulae are contained in ref. In this approach 
the deuteron form factors are expressed in terms of experimental neutron-proton (n — p) 
phase shifts in the triplet scattering channel and experimental values of nucleon EM form 
factors. 

According to ref. Q the formula for GcdiQ"^) appears as 
GUQ') = (pS20 + B22)2G«',(Q2)„ 

+ip&' + B"YGl',m- (1) 

In eq.(0) p is the constant which describes mixing of two n—p states with different orbital 
moments {I — and ^ = 2) at the point of the bound state, i.e., the deuteron. This constant 
is defined by the correspondence principle. Analysing the nonrelativistic limits of eqs.(l),(2), 
we can prove that p appears to be the standard asymptotic D/ S -ratio of the radial deuteron 
wave functions, so p = 0.0277 (numerical calculations show that the dependence of DCFF 
on the variation of p is very weak). All four elements of the matrix (s) {1,1' = 0, 2) |^ 
are taken at the bound state point s = {Md = 2M — e, where Md,M are deuteron 
and nucleon masses and e is the deuteron binding energy). All relativistic aspects of the 
two~nucleon problem are contained in G" - matrix: 



^2 r dsAB^s) n'^''*^ ds'g,{s,s't)AB{s') 

±^Vi-t)m^-tM^-^M^)- (2) 

In eq.(^ is the normalization constant, which is calculated from the condition 
Gcd{0) = 1. Matrix functions A_B" (s) = i?" (s + ie) — i?" (s — ie) are the discontinu- 
ities of the Jost matrix B{s). As usual, the Jost matrix is the solution of the boundary 
problem in two-channel scattering theory: 

^ For the choice of kinematic variables here and in eq. (2) see Appendix A. 
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Figure 1: Neutron-proton phase shifts '^Si,^ Di,^ e\ used in the calculations. Experimental 
data are taken from the VPI analysis ref.|0|. 
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The reduced Jost matrix B in eq. ( |l]) is the solution of the same eq. (||) with the scattering 
matrix S = S{5, e, ff). Expressions for B and B in terms of n—p phase shifts are cumbersome 
and are summarized in Appendix B. 

The matrix functions g]^ (s, s', t) of three variables are the relativistic charge form factors 
of the unconnected part of the matrix element of EM current {n! p' \j ^\np) . The results of 
the calculations of 5" are given in Appendix C. It is interesting to note that in the general 
case in the relativistic regime (7" - functions are not factorizable in s, s' variables, whereas 
in the nonrelativistic limit such factorization takes place. It means that in the framework 
of the used relativistic approach [Q-|6| it is impossible to introduce a concept of relativistic 
deuteron wave function. 

The experimental set of n — p phase shifts were taken from the analysis of Virginia 
Tech group and is shown in Fig. |l|. This analysis was made in the energy range 
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< Elab < 1100 McV. Extrapolation to higher energies is not as important for the calcula- 
tions of Gcd for the small and intermediate values of Q^. The only essential circumstance 
is that Si -phase shifts change sign from positive to negative and have the minimum near 
the energy Eiab ~ IGeV, then go to zero in accordance with the Levinson's theorem. Any 
realistic n — p '^Si - phase shift analysis has such a behavior. Two other states ('^-Di and 
•^Ei ) give a relatively small contribution to Gcd- 

For the calculations of Gcd we used (as a first step) the simplest choice of the nucleon 
form factors: Gep = (1 + Q^/18.23 fm~'^)''^,GMp/tJ.p = GMnl\J-n = Gsp, Geu = for aU 
Q\ 

The result of the calculations are presented in Fig. ^. Our brief conclusions are the 
following. 

1. The appearance of zero and secondary maximum in \Gcd{Q'^)\ at intermediate values 
of is the direct consequence of the change of sign of the experimental ^Si- phase shifts 
at intermediate energies. It is easy to calculate that the model's 5{E), which decreases 
monotonically with E and is always positive {5{E) > for all E), immediately leads to 
monotonically decreasing with values of \Gc'd{Q'^)\ without any fine structure. 

2. Almost up to the point of zero 

(Q2 ~ 21 fm-'^) of \Gcd{Q'^)\ the total relativistic 
correction (TRC), i.e., the difference between Gcd calculated relativistically (1,2) and its 
nonrelativistic limit, is positive and appears to be not small. For example, for ~ 14 fm~^ 
it reaches the value of 25%. 

3. TRC becomes large in the region of the secondary maximum of \Gcd{Q'^)\, increasing the 
magnitude of the form factor. 

4. The obtained results are consistent with the available data on Gcd from MIT-Bates 
Forthcoming data from Jefferson Lab E-94-018 0] are extremely important to test the 

proposed relativistic approach in the region of higher transferred momenta, where relativistic 
corrections appear to be significant. 

We would like to make the following comments to the obtained results. First, the de- 
pendence of \Gcd{Q'^)\ structure on the choice of different sets of experimental n — p phase 
shifts available from the literature is strong enough. Possible variation of 6,6,7] may shift 
the position of zero in \Gcd{Q'^)\ from the indicated point = 21 fm~'^ to the point 

— 16 fm^^ or to the point = 23 fm^^. At the same time the secondary maximum 
is located in the interval 26 < < 32 fm~^, and its height may change by a factor of 
seven. We can see that for improving our understanding of \Gcd{Q^)\ it would be desirable 
to obtain a more definite phase shifts analysis oi n — p scattering in triplet channel in inter- 
mediate energy region Eiab < 1 GeV. Secondly, let us indicate the dependence of IGcdiQ"^)] 
on the possible choice of nucleon EM form factors. Since the uncertainties of GspiQ^) in 
the considered range of are very small, the main effect in \Gcd{Q^)\ may be caused 
only by variation of GsniQ^)- It seems to be generally accepted that the maximal devi- 
ation of GEniQ^) from the zero- value approximation G£;„ = is given by known formula 
GEniQ^) = —^-nTGEpiQ^), whcre /i„ = —1.91 is the neutron anomalous magnetic moment 
and r = Q^/AAP. The results of the calculations of \Gcd{Q'^)\ with this nonzero values of 
GEniQ^) are shown in Fig. 2. One can see that the effect is sizable and the contributions of 
relativistic effects and nonzero Geu have a similar behaviour. 

Finally, we show for comparison in Fig. 2 the results of calculation of Gcd in a relativistic 
approach, developed in ref. It may be seen that zero of \Gcd{Q'^) \ predicted in ref. is 
shifted to the lower values of and the height of the secondary maximum is approximately 
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Figure 2: Relativistic deuteron charge form factor (solid line) and its nonrelativistic limit 
(dash-dotted line). A result with nonzero values of Geu = ^IJ-utGep is also shown with a 
short-dash line. A representative result of the relativistic approach of Arnold, Carlson, Gross 
1^ (dash-double-dotted line) is presented for comparison. 



the same as in our calculations. Note that in more recent calculations in the similar approach 
1^ the predicted position of zero in \Gcd{Q'^)\ remains almost unchanged. 

Here we restricted ourselves only to the discussion of the deuteron charge form factor 
Gcd ■ Even in this case we omitted such interesting questions as an analytical representation 
of relativistic corrections in different orders in (u/c)^, the new representation for realistic 
deuteron wave functions, the role of relativistic rotation of nucleon spins and orbital mo- 
mentum / = 2 in the deuteron, the problem of extraction, using the present approach, of 
GsniQ^) for ultralow values of from experimental data on elastic ed-scattering, and con- 
tributions from meson-exchange currents. It would also be interesting to perform a detailed 
comparison of the present approach with other relativistic approaches to the description of 
deuteron structure. 

All these questions, as well as the calculations of the deuteron magnetic and quadrupole 
form factors will be discussed in forthcoming publications. 
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A Kinematic variables. 

By definition s is the invariant mass oi n — p system squared: 

s = {Pn+Pp)l ■ 

In laboratory (LS) and center-of-mass (CMS) systems we have 

where E is the nucleon's energy in LS and p is modulus of the nucleon 3~momentum in 
CMS. 

is the magnitude of the 4-momcntum transfer squared: 
B Jost matrices B^B. 

The formulae for pairs {S, B) and (S, B) have the most convenient form in the p-plane: 

S{p)B^{p) = B^{p) , — CG < p < CO , 
where S = S[d{p),ri{p), e{p)], see eq.lQ). Let us introduce two new matrices S and B: 

B±{p) ^R{Tp)B±{p), 

Now the equation for B has the form 

f S{p)I3+{p) = B-{p), -oo <p<(X) , , . 

\ R{p)S{p)R-H-p)^S[S,f,,e]. ^' 

The last equation defines the reduced phase shifts S, s, fj as functions of input experi- 
mental phase shifts (5, e, ry. 

The solution of eq. (||) was found in rcf . in the form of series 

oo 

B±{p) = B±^o{p) -[I+Y. B±,'n{p)], 

m— 1 

where 

(pi(p)eT^(p) "\ 

(p2{p)e^'^''P^ J ' 

^i(p)=cxp[ V.P. — ], 

TT P' -P 



B±,o{p) = 



5 



ip2{p) = cxp[-iy.P. j 



p - p 



B±.n{p) = / ^ ■ E G,.{p')B+^o{p')[B+^,n-n{pT~'-- (6) 



In eq.(g) for odd n 
and for even n 

Gr.{p) = ^i-l)^■^mp)r■^f 0,,, 

Smn is the Kroneker delta. 



C g^} —matrix. 

In terms of invariant variables s, s' ,t and the nucleon EM form factors the matrix elements 
have the form: 

I 

g{s, s', t)[gi{s, s' , i)(cosai cosa2 — -r sinai sin 02) • GsNiQ^) + 

O 

+^^52(3, s',i) • (i sinai cosa2 - cosaisina2) • G^/jv(Q^)] : 

g{s, s',t){gi{s, s' ,t){-V2P2o cosai sina2 + -i=-P2i sinai cos 0:2) • G^^ 

V 2 

-^^52(5, s',t){V2P2o sinai cosa2 + -^-^21 cosai sina2)GXf^}, 



22/ / j.\ 

g{s, s', t){gi{s, s', t) [(^^21^21 + ^^20^20) cos(ai - as) + 

+ (^^22^22 + ■^^20^2o)cOSaiCOSa2 + 

+ (^(P22P:^i - P21P22) + ^(^2iP^o - ^20^21)) sin(ai - a2) 
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-^(-P22-P20 +-P2of'22)sinaisina2 • G'^jy- - l^92{s,s',t) 

Y^((-P21-P22 - -P22-P21) + ■^{P20P21 - ^21^20)) • 

cos(ai - 02) - (^-P22-P22 + ^-P20-P2o) COS «! sin a2 - 



-g(-P22-P20 - ^20-P22) sin ai cosa2 + 
-(^^21^21 + sin(ai - a2) 



■G 



MN 



where 



~ y/(s-4M^)(s'-4M^) ■ [A(«,s' ■ V^' 

(/i(s, s', t) = s + s' — t, 

r 1 1/2 

ff2(s,s',i)= [(-l)(M2A(s,s',t) + ss't)J , 

A(s, s', i) = + s'2 + - 2{ss' + st + s't). 
Pim are the Legendre polynomials, Pim = Pim{x) and P/^^ = P;to(x'), where 

x(s, s' , t) 



y/{s'-4:M'^)\(s,s',t) ' 

.'j;'(.s,.s',f) = -x(s',s,t). 
The angles ai,a2 of the relativistic rotation of nuclcon spins in deuteron are 



ai = arctan 
a2 = arctan 



M -t] + VF+2M) ' 

g2{s,s' ,t){y^+W+2M) 

M(s+s'-t)(Vs+vV+2M)+vW(4M2-t) ' 



T = Q^/AM"^ ; G% j^p^ = \{Ge,mp + GE,Mn) are the nucleon isoscalar charge and magnetic 
form factors. 
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